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I. INTRODUCTION 

Many measurements, such as jet quenching, elliptic 
flow, and valence quark number scaling of elliptic flow 
have provided strong evidences for the creation of the 
quark-gluon plasma (QGP) in heavy ion collisions at the 
Relativistic Heavy-Ion Collider (RHIC) and the Large 
Hadron Collider (LHC) [1- ]. With the formation of 
QGP and the restoration of chiral symmetry, strange 
and anti-strange quarks become abundant in the bulk 
medium above T c , which subsequently enhance the pro¬ 
ductions of strange and multi-strange hadrons in rela¬ 
tivistic heavy ion collisions [5]. In the past decades, dif¬ 
ferent aspects of strange and multi-strange hadrons have 
been studied in theory [5-14] and in experiment [15— 
21] . It is generally believed that multi-strange hadrons, 
such as S and LI, directly carry the information of the 
QGP phase because of their small hadronic cross sec¬ 
tions and the associated early decouplings from the sys¬ 
tem near T c [6]. Compared with common hadrons, their 
anisotropy flow are mainly developed in the QGP stage 
and less contaminated by the hadronic evolution. 

Since the running of 2.76 A TeV Pb+Pb collisions 
at the LHC, the flow and other soft hadron data of 
all charged and identified hadrons have been studied 
by many groups within the framework of hydrodynam¬ 
ics [22-31]. Using the VISHNU hybrid model [32] that 
connects (2+l)-d viscous hydrodynamics with a hadronic 
afterburner, we extracted the specific QGP shear viscos¬ 
ity ( v/ s )qgp from the elliptic data of all charged hadrons 
with MC-KLN initial conditions [24] . With the extracted 
value of (??/s)qgp> VISHNU provides a good description 
of the soft hadron data for 7r, K and p at the LHC [25]. 
Recently, the multiplicity, px-spectra and elliptic flow for 
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A, S and LI have been measured by the ALICE Collab¬ 
oration [19-21]. It is thus the right time to systemati¬ 
cally study these strange and multi-strange hadrons at 
the LHC via the VISHNU hybrid model. 

This paper is organized as follows. In Sec. II, we briefly 
introduce the VISHNU hybrid model and its setup for the 
calculations. Sec. Ill compares our VISHNU results with 
the ALICE measurements in 2.76 A TeV Pb+Pb colli¬ 
sions, including the centrality dependence of the multi¬ 
plicity density, px-spectra and differential elliptic flow for 
A, S and H. Sec. IV studies and discusses mass ordering 
of elliptic flow among n, K, p, A, 5 and Ll at the LHC. 
Sec. V explores chemical and thermal freeze-out of vari¬ 
ous hadron species during the UrQMD evolution of VISHNU. 
Sec. VI summarizes our current work and presents a brief 
outlook for the future. 


II. SETUP OF THE CALCULATION 

In this section, we describe the inputs and setup 
for the VISHNU calculations for the soft hadron data 
in 2.76 A TeV Pb+Pb collisions. The VISHNU hybrid 
model [32] combines (2+l)-d relativistic viscous hydro¬ 
dynamics (VISH2+1) [33] for the QGP fluid expansion 
with a microscopic hadronic transport model (UrQMD) [34] 
for the hadron resonance gas evolution. The transition 
from hydrodynamics to the hadron cascade occurs on 
a switching hyper-surface with a constant temperature. 
Generally, the switching temperature T sw is set to 165 
MeV which is close to the QCD phase transition tem¬ 
perature [35-37]. For the hydrodynamic evolution above 
T sw , we input an equation of state (EoS) constructed 
from recent lattice QCD data [38, 39]. 

Following Ref. [24-2 ], we input MC-KLN initial con¬ 
ditions [40, 41] and start the hydrodynamic simulations 
at tq = 0.9 fm/c. For computational efficiency, we imple- 
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ment single-shot simulations [24-26, 32, 42] using smooth 
initial entropy density profiles generated by the MC-KLN 
model through averaging over a large number of events 
within specific centrality bins 1 . Considering the conver¬ 
sion from total initial entropy to the final multiplicity of 
all charged hadrons, we cut the centrality bins through 
the distribution of total initial entropies obtained from 
the event-by-event fluctuating profiles from MC-KLN. 
Such centrality classification was once used by Shen in 
Ref. [ 4], which is more close to the experimental cut 
from the measured multiplicity distribution. The nor¬ 
malization factor for the initial entropy density is fixed 
by the charged hadron multiplicity density in the most 
central collisions ( dN c h/dr] « 1601±60 from ALICE [45]). 
The A parameter in the MC-KLN model, which quanti¬ 
fies the gluon saturation scale in the initial gluon dis¬ 
tributions [40], is tuned to 0.138 for a better fit of the 
centrality dependent multiplicity density for all charged 
hadrons. 

The QGP specific shear viscosity {t]/s)qgp is set to 
0.16 for MC-KLN initial conditions. Such combined set¬ 
ting for VISHNU once nicely described the elliptic flow 
of pions, kaons and protons in 2.76 A TeV Pb+Pb colli¬ 
sions [25]. Here, we continue to use it to further study the 
soft hadron data of strange and multi-strange hadrons at 
the LHC. To simplify the theoretical investigations, we 
set the bulk viscosity to zero and neglect the net baryon 
density and heat conductivity for the QGP systems cre¬ 
ated at the LHC. 


III. MULTIPLICITY, SPECTRA, AND 
ELLIPTIC FLOW FOR A, E AND Q 

The multiplicity, spectra and elliptic flow of pions, 
kaons and protons in 2.76 A TeV Pb+Pb collisions have 
been studied in our early paper [25]. We showed that, 
with MC-KLN initial conditions, rj/s = 0.16 and other 
parameters fixed from the related data of all charged 
hadrons, VISHNU could nicely describe the soft hadron 
data of pions, kaons and protons at the LHC. We also 
found that baryon-antibaryon (B-B) annihilations in the 
UrQMD module of VISHNU could reduce the proton yields 
by 0(30%), leading to nice fits of the proton data mea¬ 
sured by ALICE. In this section, we extend our early 
VISHNU simulations to high-statistics runs to further 
study the soft hadron data for the strange and multi- 
strange hadrons A, S and f 1 in 2.76 A TeV Pb+Pb colli¬ 
sions 2 . 


1 For recent development on event-by-event VISHNU simulations, 
please refer to [43, 44]. 

2 The 0 meson is another important multi-strange hadron that 
might directly carry the QGP information due to its small 

hadronic cross-sections. VISHNU predictions for the spectra and 
elliptic flow of 0 can be found in Ref. [25]. However, later com¬ 
parisons showed pretty large deviations between theory and ex- 



FIG. 1. (Color online) Centrality dependence of the multi¬ 
plicity density per participant pair, (dA/dy)/(A part /2) for 
A, E and Q in 2.76 A TeV Pb+Pb collisions. Insert: 
(dN/dy)/(N pa , r t / 2) for 7r, K, p and for all charged hadrons. 
Experimental data are from the ALICE Collaboration [19, 20, 
45, 46]. Theoretical curves are calculated with the VISHNU hy¬ 
brid model, using MC-KLN initial conditions, tj/s = 0.16 and 
T sw = 165 MeV. 


Figure 1 shows the centrality dependence of the mul¬ 
tiplicity density per participant pair (dN/dy)/(N palt /2) 
for A, S and Q in 2.76 A TeV Pb+Pb collisions 3 . In the 
insert, we plot the corresponding curves for pions, kaons, 
protons and for all charged hadrons that were once pre¬ 
sented in our early paper [25] for the easiness of refer- 


periment [21]. Unlike other hadrons, 0 mesons are mainly recon¬ 
structed from the strong decay channel 0 -4 K + K ~, rather than 
being directly measured. The succeeding hadronic scatterings of 
kaons might contaminate the weak signals of 0. In this paper, 
we will not show and discuss the results of 0 , but just quickly 
mention our early work [25] and leave the puzzle of 0 for future 
study. 

! We notice that the measured multiplicity of A from ALICE are 
contaminated by the feed-down decays of 1 + and E( 1385) [19]. 
However, the UrQMD module of VISHNU only includes strong reso¬ 
nance decays, but without any weak decays. To partially account 
the effects from weak decays, e.g. E° —> A + 7 , we directly sum 
the multiplicity of A and E° from VISHNU to get a corrected curve 
of A (the solid blue line with star symbols). The original yields 
of A from VISHNU are also shown in Fig. 1, which is presented by 
the dashed blue line with star symbols. In our estimations, the 
E° —> A +7 channel contributes ~ 30% additional A productions. 
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FIG. 2. (Color online) Transverse momentum spectra of A, 5 and f l at various centralities in 2.76 A TeV Pb+Pb collisions. 
Experimental data are taken from ALICE [19, 20]. Theoretical curves are calculated with VISHNU with the same inputs as for 
Fig.l. From top to bottom the curves correspond to 0-10% (xlO 4 ), 10-20% (xlO 3 ), 20-40% (xlO 2 ), 40-60% (xlO 1 ) and 60-80% 
(xl) centrality, respectively, where the factors in parentheses are the multipliers applied to the spectra for clear separation. 
Spectra of A start from 0-5% (xlO 5 ) and 5-10% (xlO 4 ), instead of 0-10%. 


ence. The inputs of our current calculations are close to 
the ones used in [25, 26], except for two points: 1) chang¬ 
ing the A parameter in the MC-KLN model from 0.28 to 
0.138, 2) cutting the centrality bins through initial en¬ 
tropy rather than the participant number N part (please 
refer to Sec. II for details). Compared with the early 
setup, these two changes mainly improve the description 
of the centrality dependent multiplicity for all charged 
hadrons and for pions, but they have small influence on 
other theoretical results, such as the elliptic flow of all 
charged and identified hadrons, etc.. 

One finds that, VISHNU nicely describes these 
(dN/dy)/(N pai t/2) curves for all investigated hadrons. 
Like the case of protons, B-B annihilations also reduce 
the yields of strange and multi-strange baryons with 
0(30%) for A, 0(20%) for 2 and f l in the most cen¬ 
tral Pb+Pb collisions(please refer to Fig. 5 in sec. V). 
The lower panel of Fig. 1 shows the difference between 
the theoretical calculated and the experimental measured 
particle yields. From the most central to semi-peripheral 
collisions, the deviations are all within 20%. For the 60- 
80% centrality bin, the differences increase to 40% for 
2, and 100% for This indicates that the strangeness 
no longer reach chemical equilibrium in the small system 
created in peripheral Pb+Pb collisions. 

In Fig. 2, we compare the transverse momentum spec¬ 
tra of A, 2 and f! from VISHNU with the measurements 
from ALICE. In general, VISHNU describes the pt spec¬ 


tra of these strange and multi-strange hadrons from the 
most central to semi-peripheral collisions, except for the 
60-80% centrality bin. Here, the theoretical curves of A 
are plotted with the original values from VISHNU with¬ 
out weak decays. As a result, they are about 30% lower 
than the experimental measurements with weak decay 
contaminations. The f l spectra from VISHNU are slightly 
higher than the experimental data for most centralities, 
but obviously above the data at the 60-80% centrality 
bin. Such deviations between theory and experiment 
are consistent with the model and data differences for 
the centrality dependent multiplicity shown in Fig. 1. In 
spite of the normalization issues, VISHNU nicely fits the 
slope of the pt spectra for A, 2 and f l at various central¬ 
ities. Together with the early nice descriptions of the pt 
spectra for pions, kaons and protons [25], it indicates that 
VISHNU generates a proper amount of radial flow, during 
its QGP and hadronic evolution, to push the spectra of 
various hadrons. 

Figure 3 presents the differential elliptic flow of A, 2 
and at three chosen centralities in 2.76 A TeV Pb+Pb 
collisions. The experimental data are from ALICE, which 
are measured with the scalar product method [21], The 
theoretical lines are calculated from VISHNU with MC- 
KLN initial conditions and (u/s)qgp = 0.16. Such in¬ 
puts once nicely described the elliptic flow data of pions, 
kaons and protons at the LHC [25]. In principle, the cur¬ 
rent calculations can be considered as extensions of the 
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FIG. 3. (Color online) Differential elliptic flow of strange hadrons A, multi-strange hadrons 5 and fl at 10-20%, 30-40% and 
50-60% centralities in 2.76 A TeV Pb+Pb collisions. Experimental data are from ALICE [21], theoretical curves are calculated 
from VISHNU with the same inputs as for Fig. 1 and Fig. 2. 



FIG. 4. (Color online) Differential elliptic flow of tt, K, p, A, S and at 10-20% and 40-50% centralities in 2.76 A TeV Pb+Pb 
collisions. Left panels (a) and (c) are plotted with V 2 from ALICE [ 1], right panels (b) and (d) are plotted with V 2 from 

VISHNU. 


early simulations [25]. Fig. 3 shows that the elliptic flow 
data below 2 GeV for A, S and fl are fairly described by 
VISHNU within the statistical error bars. Above 2 GeV, 
the descriptions of the elliptic flow for 3 at 50-60% and 
for Q at 30-40% and 50-60% become worse. On the other 
hand, viscous corrections probably become too large in 
that higher px region, making the hydrodynamic descrip¬ 
tion in VISHNU lost its predictive power. 


IV. MASS ORDERING OF ELLIPTIC FLOW 


Mass ordering of elliptic flow among various hadron 
species reflects the interaction between the radial and 
elliptic flow during the hadronic evolution. The radial 
flow tends to push the heavier particles at lower px to 
higher px, leading to a mass ordering of the px depen- 
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dent elliptic flow below 1.5-2 GeV that decreases with 
the increase of hadron mass. Such V 2 mass-ordering has 
been discovered in experiments at both RHIC and the 
LHC [21, 47-49], which has also been studied within the 
framework of hydrodynamics [11, 26, 50-52] and the blast 
wave model [53, 54]. 

In Fig. 4, we investigate mass ordering of elliptic flow 
among 7r, K, p, A, S and SI in 2.76 A TeV Pb+Pb col¬ 
lisions. For clear presentations, the ALICE data and 
the VISHNU results are plotted in separate panels for the 
two chosen centralities at 10-20% and 40-50%. Calcula¬ 
tions in [25] and in this paper (Fig. 3) have respectively 
showed that VISHNU generally describes v -2 (pr) for vari¬ 
ous individual hadrons over a wide range of centralities. 
However, further comparisons in Fig. 4 illustrate that 
VISHNU could not describe the v-i mass ordering among 
all hadron species. Although VISHNU nicely describes the 
mass ordering among 7r, K, p, and f2, it fails to correctly 
describe the mass ordering among p, A and 5. In con¬ 
trast, pure viscous hydrodynamics VISH2+1 has correctly 
predicted the relative mass-ordering among these investi¬ 
gated hadrons, but it has difficulties to roughly fit the i >2 
data for these heavier hadrons like p, A, 5 at the 10-20% 
centrality bin [51]. 

Compared with the elliptic flow of individual hadrons, 
the V 2 mass-splittings between different hadron species 
reveal more details for the hadronic evolution. Although 
VISHNU could improve the description of V 2 (pt ) for the 
hadron species like p, A and 5 through its microscopic 
hadronic scatterings, it slightly under-predicts the pro¬ 
ton V 2 below 2 GeV, leading to inverse V 2 mass ordering 
between p and A, and accidental overlaps of the elliptic 
flow for p and H below 1.5 GeV. An initial flow could 
enhance the radial flow in the hadronic stage, which is 
thus expected to improve the description of V 2 mass or¬ 
dering within the framework of the hybrid model. Mean¬ 
while, the UrQMD hadronic cross sections also need to be 
re-evaluated and improved. These have not been done 
currently and should be investigated in the near future. 


V. CHEMICAL AND THERMAL FREEZE-OUT 
OF VARIOUS HADRON SPECIES 

In this section, we investigate chemical and thermal 
freeze-out of various hadron species during the UrQMD evo¬ 
lution within the framework of the VISHNU hybrid model. 

During the QGP fireball evolution, a large number 
of hadrons are produced near T c , which subsequently 
undergo inelastic and elastic collisions in the hadronic 
phase. With the termination of inelastic collisions, the 
yields of each hadron species no longer change. The sys¬ 
tem is considered to reach chemical freeze-out. Thermal 
freeze-out happens later, which is associated with the end 
of elastic collisions. After that, the momentum distribu¬ 
tions of final produced hadrons are fixed. 

In the statistical model, the chemical freeze-out tem¬ 
perature T c h and the baryon chemical potential /i& are 


extracted from the particle yields of various hadrons [55- 
58] . A systematic study of the related data at top RHIC 
energy gives T ch ~ 165 MeV [55]. This temperature 
could describe the yields of many identified hadrons in 
2.76 A TeV Pb+Pb collisions, but obviously over-predicts 
the protons/antiprotons data at the LHC. A good de¬ 
scription of the p/p data requires a lower chemical freeze- 
out temperature around 150 MeV. However, such lower 
temperature breaks the early nice description of 5 and H 
yields once achieved with 165 MeV [59, 60]. 

To study the above proton puzzle from the statisti¬ 
cal model, we systematically investigated the soft hadron 
data for 7r, K, p at both RHIC and the LHC with the 
VISHNU hybrid model [25]. We found that baryon and 
anti-baryon annihilations influence the transport of pro- 
tons/antiprotons during the hadronic evolution, leading 
to a largely improved description of the p/p yields when 
compared with the case without B-B annihilations 4 . 
Meanwhile, other soft hadron data of n, K and p are 
also nicely fitted in general. This paper extends the early 
VISHNU calculations to further study strange and multi- 
strange hadrons at the LHC. Sec. Ill has showed a nice 
description of the paricle yields for A, 5 and f2, together 
with good fits of the spectra and elliptic flow for these 
hadrons. 

In our calculations, the switching temperature, that 
connects the hydrodynamic description of the QGP ex¬ 
pansion to the Boltzmann approach for the hadron reso¬ 
nance gas evolution, is set to 165 MeV at both RHIC and 
the LHC. However, this temperature can not be identi¬ 
fied as the chemical freeze-out temperature in the sta¬ 
tistical model, since B-B annihilations and other inelas¬ 
tic collisions are still frequent during the early hadronic 
evolution, which constantly change the yields of various 
hadrons. Instead, different hadronic scatterings in UrQMD 
lead to a hadron species dependent chemical freeze-out 
procedure. 

Figure 5 studies the time evolution of particle yield 
density for ir, K, p, A, S and f 1 during the UrQMD hadronic 
expansion. This investigation is still done within the 
VISHNU simulations for 2.76 A TeV Pb+Pb collisions, but 
exports the UrQMD intermediate results at different evo¬ 
lution times. For the easiness of comparison, Left panels 
(a)-(f) plot the time evolution of relative particle yield 
density: 0). Here, ^(t) and ^(0) denote 

the particle yield density at nrid-rapidity at later evolu¬ 
tion time and at the starting time, respectively. 

For the simulations without B-B annihilations, the 
yields of 5 and H almost do not change. This indicates 
that these two multi-strange baryons experience early 
chemical freeze-out near the switching hyper-surface of 
VISHNU. For other hadron species, their yields constantly 
change with the UrQMD evolution. By the end of the evo¬ 
lution, the yields of K and p respectively decrease ~ 5% 


4 Other related work could be found in Ref. [60—62], 
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t (fm/c) t (fm/c) 

FIG. 5. (Color online) Left panels (a)-(f): time evolution of 
relative particle yield density ^(t)//^/(0) for 7r, K, p, A, 
3 and 12 during the UrQMD expansion of VISHNU. Right pan¬ 
els (g) and (h): time evolution of the changing rate for the 
corresponding particle yield density. Solid/dashed lines de¬ 
note the VISHNU simulations with/without B-B annihilations 
in the most central 2.76 A TeV Pb+Pb collisions. 


and ~ 10%, and the yield of A increases ~ 40%. Mean¬ 
while, the changing rates for the particle yield density of 
K, p and A show wide peaks along the time axis (Panel 
(g)), illustrating that the associated inelastic collisions 
are still frequent after 10-20 fm/c. This indicates these 
hadrons experience later chemical freeze-out. The yield 
of pions only slightly increases during the UrQMD evolu¬ 
tion without B-B annihilations. However, this is not 
necessarily associated with early chemical freeze-out of 
pions. Instead, pions maintain relative chemical equilib¬ 
rium below T c through frequent quasi-elastic collisions, 
e.g., 7T7r o p, 7 tN o A, etc. [63-66]. 

The B-B annihilations (pp —> mr, etc.) mainly influ¬ 
ence the baryon’s transport in UrQMD, leading to ~ 30% 
reductions for the p and A yields, ~ 20% reductions for 
the S and 12 yields in the most central Pb+Pb collisions. 
Meanwhile, the yields of 7r and K slightly increase by 
~ 5% through the annihilation channels. Although these 
two multi-strange hadrons, S and 12, rarely interact with 
other hadrons during the UrQMD evolution, the annihila¬ 
tions with their own anti-particles delay their chemical 
freeze-out. This is presented by the wide peaks on the 
changing rate curves for these two multi-strange baryons 



t (fm/c) t (fm/c) 

FIG. 6. (Color online) Similar to Fig. 5, but for 70-80% cen¬ 
trality bin. 


in panel (h). But, compared with other curves such as 
the proton one, E and 12 still experience early chemical 
freeze-out. The B-B annihilations almost balance with 
the other inelastic collision channels on the production 
of A and K. As a result, the yields of these two hadrons 
only slightly change during the hadronic evolution. 

Fig. 6 is similar to Fig. 5, but for the 70-80% centrality 
bin. For the case without B-B annihilations, the parti¬ 
cle yields of various hadrons almost do not change dur¬ 
ing the UrQMD evolution. Compared with the most cen¬ 
tral Pb+Pb collisions, the number of in-elastic collisions 
in the hadronic phase are greatly reduced. Fig. 6 also 
showed that the B-B annihilations decrease the baryon 
yields for p, A, S and 12 by 5-8% in peripheral collisions, 
but most of the annihilations happen before 10 fm/c. 

Shortly speaking, Fig. 5 and Fig. 6 mainly concen¬ 
trate on studying time evolution of various hadron yields, 
which indirectly reflect the inelastic collisions in UrQMD. A 
further analysis of the space-time distributions of the last 
inelastic collisions will reveal direct information for chem¬ 
ical freeze-out, which may even help us to extract effec¬ 
tive chemical freeze-out temperatures of various hadron 
species. Unfortunately, the current version of UrQMD does 
not record such intermediate information. We have to 
leave it for future study. 

Besides 4-momentum of final produced hadrons, UrQMD 
also outputs the positions (in space and time) of the last 
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FIG. 7. (Color online) Thermal freeze-out time distributions for 7r, K, p, A, E and 1 1 in the most central Pb+Pb collisions, 
calculated from hydrodynamics (green), hydrodynamics+resonance decay (blue), VISNU without B-B annihilations (red) and 
with B-B annihilations (black). 


elastic collisions or resonance decays that directly reflect 
thermal freeze-out of the evolving system. Here we de¬ 
fine the time distributions of the last collisions or decays 
for various hadrons species as the corresponding thermal 
freeze-out time distributions. On the other hand, they 
can also be considered as the production-time distribu¬ 
tions of specific hadron species during the UrQMD evolu¬ 
tion. 

Figure 7 shows thermal freeze-out time distributions 
for 7r, K, p, A, 3 and SI in the most central Pb+Pb col¬ 
lisions. To study the hadronic scattering effects, we set 
four comparison simulations: 1) viscous hydrodynam¬ 
ics terminated at T sw = 165 MeV with only thermal 
hadron emissions, 2) viscous hydrodynamics terminated 
at T sw with thermal hadron productions and succeeding 
resonance decays, 3) VISHNU without B-B annihilations, 
4) VISHNU with B-B annihilations. Here, both viscous 
hydrodynamics and VISHNU simulations input the same 
initial conditions, EoS, and other related parameters as 
described in Sec. II. 

The thermal freeze-out time distributions for various 
hadrons in case (1) all stop around 10 fm/c, because 
the hydrodynamic evolution terminates around that time 
in the most central collisions. Comparing the thermal 
freeze-out time distributions from hydrodynamics (case 
1) and from hydrodynamics+resonance decays (case 2), 
we find a certain portion of the resonance decays happen 
near the hydrodynamic freeze-out surface, which largely 
enhance the productions of pions and protons before 10 
fm/c. Meanwhile, the long-lived resonances also con¬ 
tribute later hadron productions after 10 fm/c, which 


results in long tails for the distribution curves of 7r, K, 
p, A and 3. We notice that there is no change for the 
17 curves between case (1) and case (2). UrQMD only in¬ 
cludes hadrons below 2 GeV, the associated resonance 
decays do not contribute to the production of this heavy 
multi-strange baryon. 

The UrQMD hadronic scatterings in case (3) and (4) 
broaden thermal freeze-out time distributions of all 
hadron species, which shift the averaged hadron produc¬ 
tion times before 10 fm/c in case (1) and (2) to later 
values ranging from 10 fm/c to 40 fm/c. We also observe 
that the B-B annihilations further decrease the produc¬ 
tions of p, A, 3 and 17 as shown in Fig. 5. In general, such 
annihilations do not change the shape of these thermal 
freeze-out time curves. 

Figure 8 compares the thermal freeze-out time distri¬ 
butions from the most central collisions with the ones 
from peripheral collisions. Here, the results are from the 
VISHNU simulations with B-B annihilations. We find that 
the peaks of the thermal freeze-out time distributions for 
7r, K, p, A, 3 and 17 are all shifted to much earlier time 
in peripheral collisions, because the created QGP fireball 
there has much smaller volume and shorter lifetime. 

Figure 9 compares the thermal freeze-out time dis¬ 
tributions for 7T, K, p, A, 3 and 17 in the most cen¬ 
tral Pb+Pb collisions. These curves are calculated from 
VISHNU with B-B annihilations, which are the same as 
the corresponding ones shown in different panels of Fig. 7. 
The peaks of the 3 and 17 curves are both located around 
10 fm/c. Compared with these curves of p and A, which 
peaks are located around 20-30 fm/c, these two multi- 
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FIG. 8. (Color online) Comparisons of the thermal freeze-out time distributions for n, K, p, A, S and Cl between the most 
central and peripheral Pb+Pb collisions, calculated from VISHNU with B-B annihilations. 
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FIG. 9. (Color online) Thermal freeze-out time distributions 
for 7r, K, p, A, 3 and in the most central Pb+Pb collisions, 
calculated from VISHNU with B-B annihilations. 


strange hadrons experience earlier thermal freeze-out. 
Although the peaks of the n and K curves are closer 
to the ones of 5 and Cl, their freeze-out time distribu¬ 
tions spread widely along the time axis. This indicates 


that these two meson species still surfer a certain amount 
of hadronic scatterings even during the late evolution 
of UrQMD. We conclude from Fig. 9 that thermal freeze- 
out is hadron species dependent. Compared with other 
hadrons, the two multi-strange hadrons 3 and Cl experi¬ 
ence earlier thermal freeze-out, as expected, due to their 
much smaller hadronic cross sections. 


VI. SUMMARY AND OUTLOOK 

In this paper, we studied the soft hadron data of 
strange and multi-strange hadrons at the LHC, using 
VISHNU hybrid model. We found that, with MC-KLN ini¬ 
tial conditions, 77 /s = 0.16 and other inputs that fit the 
related data of common hadrons [25], VISHNU generally 
describes the multiplicity, px-spectra and differential el¬ 
liptic flow of the strange hadron A and the multi-strange 
hadrons 5 and SI in 2.76 A TeV Pb+Pb collisions. 

Compared with the pure hydrodynamic calculations 
from VISH2+1 [11], VISHNU improves the descriptions of 
the elliptic flow for p, A, 5 and Cl with microscopic 
hadronic scatterings that re-balance the interactions be¬ 
tween radial and elliptic flow. However, mass ordering 
of elliptic flow among various hadron species is not fully 
described. VISHNU slightly under-predicts the differential 
elliptic flow of protons, leading to inverse mass-ordering 
among p, A and 5. An initial flow and improved hadronic 
cross sections in UrQMD may help to solve this issue. This 
should be investigated in the near future. 

With a nice description of the particle yields for 7 r, K, 
p, A, 5 and Cl, we further investigated chemical and ther- 
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mal freeze-out of various hadron species within the frame¬ 
work of VISHNU hybrid model. We found that, compared 
with other hadrons, the two multi-strange hadrons, H and 

experience earlier chemical and thermal freeze-out due 
to their small hadronic cross sections. A study for time 
evolution of the hadron yields also shows that baryon- 
antibaryon annihilations in UrQMD delay the chemical 
freeze-out of S and fl, when compared with the case with¬ 
out such annihilations. 

We also emphasized that the switching temperature in 
VISHNU could not be identified as the chemical freeze- 
out temperature in the statistical model since inelastic 
collisions are still frequent during the early evolution of 
UrQMD. A further analysis of the space-time distributions 
of the last inelastic collisions could reveal more informa¬ 
tion on chemical freeze-out, which may even help us to 


extract effective chemical freeze-out temperatures of var¬ 
ious hadron species. Unfortunately, such investigation 
could not be done with the current version of UrQMD. With 
an updating of UrQMD to further record the intermediate 
evolution information, the chemical freeze-out procedure 
of the evolving hadronic system will be further studied 
in the future. 
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